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ABSTRACT 


The  arrivals  of  P-waves  at  the  Large  Aperture  Seismic  Array 
(LASA)  from  over  3000  earthquakes  at  various  azimuths  and  dis¬ 
tances  have  been  studied.  From  each  set  of  arrival  times  a 
variety  of  parameters  related  to  the  travel  time  curve  are 
determined.  This  report  consists  primarily  of  a  summary  of  the 
principal  features  of  these  parameters,  as  functions  of  epi- 
central  distance  and  azimuth. 
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INTRODUCTION 


Investigations  ol  the  arrival  times  of  P-waves  at  the  Large 
Aperture  Seismic  Array  (LASA)  at  Billings,  Montana,  have  been 
carried  out  by  the  principal  investigator  during  the  period 
October,  1965  to  July,  1971,  under  several  ARPA  contracts. 

These  studies  were  unusual  in  that  a  very  large  quantity  of 
high-quality  data  were  obtained.  Portions  of  these  data,  together 
with  interpretations,  have  been  published  (see  References),  but 
inevitably  much  remains  that  cannot  be  interpreted  at  present. 

The  purpose  of  this  report  is  to  catalogue  as  many  as  possible 
of  the  observations  that  have  been  made. 

Since  many  similar  studies  using  various  arrays  have 
appeared,  and  continue  to  appear,  in  the  literature,  it  is 
worthwhile  to  emphasize  the  uniqueness  of  the  present  set  of 
data.  The  quality  of  the  data  is  primarily  due  to  the  physical 
size  (about  200  km  aperture)  of  the  LASA  array,  and  to  the 
availability  of  recordings  that  permit  the  relative  timing  of 
arrivals  to  .05  seconds  or  better.  This  immediately  suggests 
that  it  should  be  possible  to  determine  dT/dA  from  these  arrivals 
with  a  precision  of  better  than  0.1  seconds/degree,  and  the 
observed  scatter  in  the  determinations  agrees  wit!  this  conclusion. 
In  fact,  the  observed  scatter  is  between  3  and  10  times  smaller 
than  those  of  nost  other  published  studies. 

The  high  quantity  of  data  is  partly  due  to  the  sensitivity 
of  the  array,  and  partly  due  to  its  location.  The  array  lies 
in  such  a  position  that  most  of  the  circum-Pacif ic  belt  of 
seismic  activity  lies  within  two  rather  small  azimuth  ranges 
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from  the  array  (see  Section  2) .  The  array  is  therefore  partic¬ 
ularly  suited  to  observations  of  the  variation  of  travel-time 
parameters  with  epicentral  distance.  It  is  unfortunate  that 
the  azimuthal  distribution  of  events  is  much  less  complete.  One 
important  result  of  the  large  number  of  earthquakes  studied  is 
that  it  has  become  possible  to  apply  strong  criteria  in  the 
selection  of  data.  The  data  shown  in  this  report  represent  only 
about  70%  of  those  obtained;  the  rest  were  eliminated  as  being 
somewhat  less  reliable  for  a  variety  of  reasons. 

One  unexpected  problem  that  results  from  the  high  quality 
of  the  present  data  is  that  considerable  fine  structure  becomes 
apparent.  This  is  shown  not  only  by  the  distance  variation  of 
dT/dA ,  but  also  in  the  differences  between  observations  along 
different  azimuths.  Unfortunately  it  is  not  possible  to  provide 
a  convincing  explanation  for  most  of  these  features.  On  the 
other  hand,  smoothing  (or  averaging)  of  the  data  does  not  seem 
justified  in  view  of  the  small  amount  of  scatter.  This  leads 
to  certain  ambiguities  in  the  interpretation  of  the  data. 

This  report  concentrates,  therefore,  on  an  exposition  of 
the  data,  and  the  methods  by  which  they  were  obtained.  Very 
few  comments  on  the  possible  meanings  of  the  results  are  included. 
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I.  ANALYSIS  OF  THE  DATA 

The  data  used  in  this  study  consists  of  the  first  arrivals 
of  P-waves  at  the  Large  Aperture  Seismic  Array  (LASA)  in  Billings, 
Montana.  These  arrivals  were  determined  visually  from  the  16  mm 
film  records  of  the  short  period  instruments  at  the  array.  In 
particular,  the  arrivals  were  recorded  for  the  center  seismometer 
in  each  of  the  sub-arrays  in  the  E  and  F  rings ,  and  the  center 
sub-array  AO.  The  locations  of  these  sub-arrays  are  shown  in 
Figure  1,  and  their  coordinates  with  respect  to  a  N-S  E-W  coordinate 
system  centered  on  the  sub-array  AO  are  shown  in  Table  1. 

The  recording  of  a  typical  earthquake  (one  of  magnitude  4.4 
in  the  Aleutian  Islands)  is  shown  in  Figure  2a.  The  lack  of 
coherence  of  the  waveforms  beyond  the  first  half-cycle  is  clearly 
indicated.  For  this  reason,  timing  methods  that  involved  cor¬ 
relation  Af  the  waveform  across  the  array  were  tested  and  eliminated 
early  in  the  study.  Ideally,  the  "knee"  of  the  onset  should  be 
timed,  but  this  was  found  to  be  difficult  for  all  but  the  largest 
events.  As  a  compromise,  all  arrival  measurements  were  made  on 
the  first  pe<  k  (or  trough)  of  the  arrival.  For  all  the  events 
included  in  this  report,  it  was  possible  to  time  this  point  with 
a  precision  of  0.05  seconds  or  better.  The  relative  times  of 
arrival  with  respect  to  AO  were  therefore  subject  to  a  maximum 
error  of  ±0.05  seconds,  and  in  most  cases  the  estimated  error  is 
about  one  half  of  this. 

Many  events  were  not  recorded  at  all  nine  seismometers,  due 
to  equipment  failure,  excessive  noise  level,  and  so  on.  It  was 
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Fig.  1:  Location  and  deonetry  of  the  largo  Aperture 
Jeisric  array  (Ln.Li) ,  near  dillingn,  on  tana. 
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found  in  the  course  of  the  study  that  the  dT/dA  values  for  these 
events  showed  a  significantly  larger  scatter  than  for  those  recorded 
at  all  the  sites.  For  this  reason,  incompletely  recorded  events 
have  not  been  included  in  this  study. 

The  basic  data  for  each  event  therefore  consists  of  a  set  of 
eight  times,  relative  to  AO.  These  times  were  processed  on  the 
IBM  360  computer  at  Brown  University  to  extract  a  variety  of  infor¬ 
mation.  One  page  of  printed  output  was  obtained  for  each  event, 
and  a  typical  output  is  shown  in  Figure  2b,  for  the  event  shown  in 
Figure  2a.  The  various  items  on  this  output  are  described  briefly 
below: 

a)  Using  the  USCGS  epicenter,  the  azimuth  and  distance  of  the 
event  were  calculated,  using  a  sub-program  that  includes  the 
effect  of  ellipticity. 

b)  The  distance  was  corrected  for  the  effect  of  focal  depth,  by 
extrapolating  the  ray  path  back  to  the  earth's  surface.  This 
correction  is  small  for  events  at  teleseismic  distances  with 
focal  depths  less  than  about  200  km.  It  can  be  estimated 
accurately  enough  for  the  present  purposes  on  the  basis  of  a 
standard  crust  and  upper  mantle  structure.  The  Jeffreys- 
Bullen  structure  was  used,  and  the  corrections  for  various 
focal  depths  and  distances  are  given  in  Figure  3.  Earthquakes 
deeper  than  200  km  were  not  used  in  this  study. 

c)  The  dT/dA  corresponding  to  any  given  model  is  determined  for 
the  corrected  distance,  and  listed  as  "calculated  slope."  In 


this  case,  the  model  A  (Chinnery,  1969)  was  used  as  a  standard. 
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This  facilitates  the  comparison  of  this  particular  event  with 
other  results. 

d)  A  plane  wave  is  first  fitted  to  the  wavefront.  If  the  arrival 
time  (relative  to  AO)  is  at  a  site  with  coordinates  x^, 
(again,  relative  to  AO) ,  the  data  are  fitted  to  an  equation  of 
the  form: 

T  +<'<  ■*-.  +  k  ;-jt  r  O 

1  L 

(1) 


The  coefficients  a  and  b  were  determined  oy  a  least  squares 
process,  and  from  them  the  apparent  azimuth  and  a  first  estimate 
of  dT/dA  could  be  obtained.  These  results  are  listed  under  the 
heading  "plane  wave  fit."  The  agreement  between  the  apparent 
and  true  azimuths  was  in  general  very  good.  This  point  is  dis¬ 
cussed  further  later  in  this  report, 
e)  Because  of  the  comparatively  large  aperture  of  the  array,  about 
200  km,  it  was  necessary  to  include  the  possibility  of  curvature 
of  the  wavefront  in  the  vertical  plane  containing  the  azimuth  to 
the  source.  A  simple  and  effective  way  of  doing  this  is  to 
project  the  positions  of  the  sub-arrays  onto  this  plane,  and 
replace  the  coordinates  x^,  y^  by  a  single  variable  the 

distance  of  the  ith  sub-array  from  AO  along  the  azimuth  towards 
the  event.  In  this  calculation,  the  "true"  azimuth  was  used. 

Now  T^  may  be  regarded  as  a  function  of  alone,  and  first, 

second  and  third  order  relations  may  be  written  as: 


r 

-r-  A  o  L 


I 


o 
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T.  -  T„  +  b  S, 

~i~,  -  T„  r  d  Sc 

where  a-f  are  constant  coefficients,  and  Tq  is  included  so 
that  the  equations  are  not  constrained  to  pass  through  AO. 

Equations  2,  3  and  4  were  fitted  to  the  observed  times 
by  least  squares,  and  it  is  clear  that  the  first,  second  and 
third  derivatives  of  the  travel  time  curve  may  be  determined 
from  the  coefficients  a-f.  The  results  are  listed  under 
first,  second  and  third  degree  polynomial  fit,  on  the  output 
sheet.  In  addition,  an  equation  of  the  form  of  4,  but  with 
d  constrained  to  be  equal  to  b,  was  used.  The  latter  fit 
was  found  to  give  the  least  scatter  in  estimates  of  the 
second  derivative. 

f)  At  the  foot  of  the  page,  the  observed  times,  the  &  ^ ,  and  the 
residuals  from  the  various  fits  are  listed.  In  addition,  the 
times  predicted  by  the  model  and  the  apparent  station  cor¬ 
rection  are  included. 

g)  Also  printed  out  are  the  differences  between  the  true  and 
observed  azimuths,  and  between  the  observed  and  calculated 
dT/d A  ,  together  with  the  standard  deviation  of  the  various 
polynomial  fits. 
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II.  DISTRIBUTION  OF  EVENTS 

Over  three  thousand  events  have  been  processed  in  the  above 
manner.  This  covers  the  operation  of  the  array  from  its  inception, 
in  October,  1965,  through  the  end  of  1968.  These  events  are  dis¬ 
tributed  at  a  variety  of  azimuths  and  distances  from  the  array, 
as  shown  in  Figure  4.  It  should  be  noted  that  only  those  earth¬ 
quakes  that  gave  an  amplitude  of  greater  than  5  millimicrons  at 
the  array  have  been  included.  This  means  that  there  is  a  con¬ 
siderable  variation  in  magnitude  of  the  events  considered,  and 
accounts  for  the  scarcity  of  data  at  distances  greater  than  about 
100°. 

It  is  clear  from  Figure  4  that  most  of  the  events  recorded 

lie  within  two  azimuth  ranges.  A  large  number  of  events  were 

recorded  from  the  northwest  of  the  array  (azimuth  range  290  to 

330) ,  and  a  slightly  smaller  number  were  recorded  from  the  south 

o  o 

(azimuth  range  140  to  170)  .  In  these  two  ranges  of  azimuth  \t 
has  been  possible  to  study  the  variation  of  dT/dA  with  distance, 
in  some  detail.  Information  is  available  at  some  other  azimuths, 
but  the  coverage  is  much  less  complete. 

It  should  be  noted,  too,  that  the  earthquakes  are  not  well 
distributed  in  azimuth.  One  of  the  interesting  results  of  this 
study  is  that  there  are  significant  azimuthal  variations  in  the 
observed  quantities.  However,  the  gaps  in  azimuthal  coverage 
make  it  difficult  to  study  these  azimuthal  variations  in  detail 
(see  section  $) . 
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III .  OBSERVATIONS  OF  dT/dA 

In  this  section  we  show  the  raw  dT/dA  data,  obtained  ai 
described  above.  The  particular  value  plotted  is  the  derivative 
resulting  from  the  second  degree  polynomial  fit  (Figure  2b) .  It 
should  be  emphasized  that  none  of  these  data  have  been  corrected 
for  source  or  array  structure.  The  nature  of  these  corrections 
will  be  considered  later. 

In  order  to  minimize  the  effects  of  the  array,  the  data 
are  presented  in  a  series  of  azimuth  ranges.  The  differences 
between  these  azimuth  ranges  are  considered  in  section  «1,  Some 
comments  are  made  on  unusual  features  in  each  azimuth  range. 

a.  Azimuth  Range  290-330  (Figures  5  and  6) 

This  azimuth  range  contains  the  most  complete  coverage  as 
a  function  of  distance.  The  scatter  of  the  data  is,  in  general, 
quite  small,  of  the  order  of  0.1  seconds/degree.  In  two  regions, 
one  between  33  and  38  degrees  and  one  between  84  and  88  degrees, 
the  scatter  is  significantly  larger  than  this.  As  we  shall  see 
in  sections  7  and  there  are  reasons  for  proposing  that  the 

first  of  these  is  the  result  of  mantle  structure,  and  the  second 
is  due  to  source  effects.  The  interpretation  of  this  azimuth 
range  is  taken  up  in  detail  in  section  7, 

b.  Azimuth  Range  330-60  (figures  7  and  8) 

The  distribution  of  events  in  this  azimuth  range  is  much 
poorer  than  above.  The  data  at  distances  beyond  75  degrees  show 
little  dependence  on  azimuth.  They  do,  however,  show  a  distinctly 
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III.  OBSERVATIONS  OF  dT/dA 

In  this  section  we  show  the  raw  dT/dA  data,  obtained  as 
described  above.  The  particular  value  plotted  is  the  derivative 
resulting  from  the  second  degree  polynomial  fit  (Figure  2b)  .  It 
should  be  emphasized  that  none  of  these  data  have  been  corrected 
for  source  or  array  structure.  The  nature  of  these  corrections 
will  be  considered  later. 

In  order  to  minimize  the  effects  of  the  array,  the  data 
are  presented  in  a  series  of  azimuth  ranges.  The  differences 
between  these  azimuth  ranges  are  considered  in  section  3.  Some 
comments  are  made  on  unusual  features  in  each  azimuth  range. 

a.  Azimuth  Range  290-330  (Figures  5  and  6) 

This  azimuth  range  contains  the  most  complete  coverage  as 
a  function  of  distance.  The  scatter  of  the  data  is,  in  general, 
quite  small,  of  the  order  of  0.1  seconds/degree.  In  two  regions, 
one  between  33  and  38  degrees  and  one  between  84  and  88  degrees, 
the  scatter  is  significantly  larger  than  this.  As  we  shall  see 
in  sections  7  and  S,  there  are  reasons  for  proposing  that  the 
first  of  these  is  the  result  of  mantle  structure,  and  the  second 
is  due  to  source  effects.  The  interpretation  of  this  azimuth 
range  is  taken  up  in  detail  in  section  7, 

b.  Azimuth  Range  330-60  (Figures  7  and  8) 

The  distribution  of  events  in  this  azimuth  range  is  much 
poorer  than  above.  The  data  at  distances  beyond  75  degrees  show 
little  dependence  on  azimuth.  They  do,  however,  show  a  distinctly 
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different  behavior  from  those  along  the  290-330  direction.  In  the 
distance  range  43  to  50  degrees,  events  from  an  azimuth  of  330-60 
depart  strongly  from  those  at  an  azimuth  of  20  to  60  degrees. 

This  may  be  the  result  of  source  structure,  but  the  data  are  suf¬ 
ficiently  sparse  that  it  is  hard  to  investigate  the  problem  further 
at  present. 

c.  Azimuth  Range  60-130  (Fig  .res  9  and  10) 

These  data  too  are  quite  sparse.  This  does  seem  to  be  a 
significant  difference  between  events  at  110-130  and  those  at 
60-110,  in  the  distance  range  45  to  60  degrees. 

d.  Azimuth  Range  130-170  (Figures  11  and  12) 

Good  coverage  is  obtained  in  this  azimuth  range.  Some  of 
the  apparent  scatter  in  this  diagram  may  be  due  to  errors  in 
epicentral  location.  The  pronounced  difference  between  the  shape 
of  this  curve  and  that  for  the  azimuth  range  290-330  is  discussed 
in  section 

e.  Azimuth  Range  170-290  (Figures  13  and  14) 

These  events  are  very  poorly  distributed  in  distance. 
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IV.  OBSERVATIONS  OF  d2T/dA2 

The  second  derivative  of  the  travel  time  curve  can  easily 
be  determined  from  the  analysis  outlined  in  section  I ,  and  is 
included  in  the  printed  output  for  each  event  (Figure  2b) .  This 
is  easiest  seen  if  the  time  T^  at  the  ith  station  is  related  to 
the  distance  of  the  ith  station  from  AO  along  the  event 
azimuth  by  a  Taylor  expansion 


d^T 

Estimates  of  g-^-2  may  therefore  be  obtained  from  the  coefficients 
d  and  f  in  equations  3  and  4 . 

This  is  the  first  time  that  an  attempt  has  been  made  to 
directly  measure  the  second  derivative  of  the  travel  time  curve 
using  an  array  (Chinnery  1968a,  1968b) .  The  results  are  shown, 
for  the  two  principal  azimuth  ranges,  in  Figures  15-18. 

It  is  clear  that  the  scatter  in  these  data  is  large,  and 
this  is  primarily  due  to  the  timing  accuracy.  The  absolute  value 
of  d  T/cfA  is  quite  unreliable,  since  station  corrections  have 
not  been  applied,  and  because  these  corrections  could  not  be 
estimated  accurately  enough  anyway  for  the  present  purpose . 
However,  there  is  one  feature  of  interest  in  these  curves,  for 
in  the  vicinity  of  a  cusp  or  triplication  of  the  travel  time  curve 
the  magnitude  of  the  second  derivative  should  become  large.  Some 
indication  of  possible  cusps  may  be  obtained  at  those  distance 
ranges  where  the  scatter  of  these  graphs  becomes  large.  These 


are  noted  below: 
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a.  Azimuth  Range  290-330  (Figures  15  and  16) 

2  2 

Unusually  large  or  small  values  of  d  T/dA  are  noted  in 

the  distance  ranges  26-27°,  36-37°  (only  two  points) ,  60-65°, 

68-69°,  75-76°,  86-90°  and  96-97°.  There  is  a  high  correlation 

between  these  anomalous  regions  and  those  selected  on  the  basis 

of  the  dT/d A  data  (section  7).  This  suggests  that  there  is  some 

2  2 

useful  information  in  the  d  T/dA  data,  though  the  present 
scatter  makes  any  interpretation  of  the  second  derivative  data 
alone  quite  meaningless. 

b.  Azimuth  Range  130-170  (Figures  17  and  18) 

Anomalous  distance  ranges  are  32-34°,  49-50°,  65-67°  and 
81-83°.  The  poor  correlation  of  these  with  the  290-330  azimuth 
range  is  striking.  We  shall  see  later  that  the  dT/dA  data  for 
these  two  azimuth  ranges  also  show  little  correlation  (section  &) . 
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V.  OBSERVATIONS  OF  AZIMUTH  ERROR 

Two  estimates  of  the  azimuth  to  an  event  may  be  obtained. 

The  "true"  azimuth  may  be  obtained  from  the  published  USCGS 
epicenter  for  each  event.  The  error  in  this  azimuth  from  its 
real  value  is  likely  to  vary  from  place  to  place,  and  may  contain 
systematic  components  which  may  create  apparent  variations  with 
azimuth  and/or  distance.  However,  the  error  in  the  "true" 
azimuth  is  likely  to  be  small.  Most  events  should  be  located 
to  within  1°  or  less  of  their  true  epicenters ,  and  the  corres¬ 
ponding  error  in  azimuth  should  be  much  less  than  this. 

A  second  estimate,  the  "observed"  azimuth,  may  be  obtained 
from  the  coefficients  a  and  b  in  equation  1.  The  difference 
between  the  true  and  observed  azimuths  is  a  measure  of  the 
bending  of  the  ray  path  out  of  th*  vertical  plane  containing  the 
source  and  the  array.  However,  where  this  bending  occurs  is 
harder  to  say.  It  may  be  due  to  source  structure,  a  lateral 
variation  in  mantle  structure,  or  local  structure  under  the  array. 

Figure  19  shows  the  variation  of  true  minus  observed  azimuth 
with  distance  for  events  to  the  Northwest  of  the  array.  Clearly 
there  is  a  strong  but  apparently  continuous  variation  of  this 
quantity.  The  continuity  suggests  that  the  variation  is  not  due 
to  source  structure,  though  this  cannot  be  proved.  The  large 
errors  for  events  at  distances  greater  them  100°  is  very  inter¬ 
esting,  since  these  rays  should  trrverse  crustal  structure  almost 
vertically,  and  they  should  therefore  be  least  affected  by  the 


crustal  structure. 
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Figure  20  shows  similar  data  from  the  Southeast  events.  The 
trend  of  the  observations  shows  a  marked  similarity  with  Figure  19 
in  the  distance  range  50-90°,  but  is  very  different  outside  this 
range . 

One  possible  reason  for  an  azimuth  error  is  crustal  structure 
beneath  the  array.  If  this  structure  may  be  approximated  by  a 
dipping  interface,  the  graph  of  azimuth  error  against  azimuth 
should  be  a  sine  function  (Niazi,  1966).  Figures  21  and  22  show 
such  graphs,  for  events  in  the  distance  ranges  20-55°  and  75-95° 
respectively.  No  clear  sinusoidal  terms  are  apparent,  and  in 
fact  there  are  significant  differences  between  the  two  sets  of 
data. 

It  has  not  proved  possible,  as  yet,  to  interpret  these 
observations  of  azimuth  error. 
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VI.  OBSERVATIONS  OF  TRAVEL-TIME  ANOMALIES 

There  is  no  doubt  that  the  relative  arrival  times  of  P-waves 
across  LASA  show  a  surprising  complexity.  Many  studies  (e.g., 
Chiburis,  1966)  have  confirmed  this,  and  there  is  little  doubt 
that  this  complexity  is,  at  least  in  part,  due  to  structure 
beneath  the  array.  However,  attempts  to  use  this  complexity  to 
determine  the  structure  have  met  with  little  success. 

At  least  part  of  the  reason  for  this  lies  in  the  problem  of 
defining  the  nature  of  the  travel  time  anomalies  at  each  sub-array. 
There  are  two  basic  ways  to  define  these  anomalies.  First,  the 
arrival  times  at  the  sub-arrays  can  be  compared  with  those  pre¬ 
dicted  by  a  standard  travel-time  curve.  Second,  the  arrivals 
themselves  may  be  fitted  with  a  straight  line  or  a  curve,  and 
the  residuals  from  the  fit  defined  as  travel  time  anomalies. 

Each  of  these  definitions  has  advantages  and  disadvantages,  and 
may  be  applied  to  the  data  in  various  ways.  They  are  discussed 
below: 

(i)  Departure  from  Standard  Travel  Times. 

If  the  distance  of  an  earthquake  from  the  array  is  known, 
then  any  standard  travel  time  table  may  be  used  to  estimate 
the  relative  arrival  times  at  the  various  sub-arrays.  Many 
authors  have  used  the  Jef freys-Bul len  Tables  for  this  purpose 
(e.g.,  Chiburis,  1966,  Lincoln  Laboratories  Report  LL-6,  1967, 
Greenfield  and  Sheppard,  1968) .  Other  possiblities  may  include 
the  Herrin  Tables  or  any  velocity  model  for  the  earth's  interior, 
such  as  the  model  A  of  Chinnery  (1969)  (see  section  7) . 


It  would  be  nice  if  we  could  define  the  travel  time  anomaly 


at  the  ith  station  by 

/W' 


(6) 


P  o 

where  T.  is  the  predicted  time  and  lb  is  the  observed  time.  We 

cannot  do  this,  of  course,  since  (because  of  origin  time  errors) 

the  absolute  value  of  T.°  is  not  known.  We  are  forced  therefore 

1 

to  select  some  reference  time  against  which  to  measure  the  T^°. 
The  simplest  way  to  do  this  is  to  choose  the  arrival  at  one  sub¬ 
array  as  a  reference,  and  write 


At,  -  ZHt. 


(7) 


The  right  hand  side  of  this  equation  may  be  determined  easily 
from  the  observations,  and  many  authors  have  defined  the  "travel 
time  anomaly"  at  the  ith  station  as  this  quantity  (Chiburis,  1966, 
Lincoln  Laboratories  Report  LL-6  ,  1967  ,  etc.).  This  definition 
may  be  satisfactory  for  calibrating  the  array,  but  is  quite 
inadequate  for  studies  of  the  crustal  structure  at  the  array, 
for  two  reasons  . 

In  the  first  place,  this  definition  includes  only  the  dif¬ 
ference  between  the  structural  contributions  at  the  ith  station 
and  the  reference  station.  This  is  particularly  important  if 
the  station  AO  is  chosen  as  reference,  since  AO  is  one  of  the 
more  complex  sites.  Secondly,  this  definition  of  travel  time 
anomaly  attributes  any  difference  between  the  standard  model  and 
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the  true  travel  time  curve  to  an  effect  at  the  sub-arrays.  There 
is  much  evidence  that  any  world-wide  average  travel  time  table 
may  depart  significantly  from  measurements  along  an  individual 
profile  on  the  earth's  surface,  and  this  is  particularly  true 
of  the  rather  out-dated  Jef f reys-Bullen  Table.  Each  "travel 
time  anomaly"  will  therefore  be  contaminated  both  by  structure 
at  the  reference  site  and  structure  in  the  earth's  interior. 

The  second  of  these  problems  can  be  minimized  by  using  a 
set  of  travel  times  determined  directly  from  the  LASA  data.  As 
an  example,  Figure  23  shows  the  travel-time  anomalies  for  sub¬ 
array  Fl  calculated  using  equation  7,  with  AO  as  a  reference  sub¬ 
array,  and  predicted  travel  time  differences  according  to  model 
A  (Chinnery,  1969;  section  7  of  this  report).  The  anomalies  for 
events  to  the  Northwest  of  the  array  are  plotted  as  a  function 
of  epicentral  distance,  and  show  a  strong  variation  with  distance. 

The  contribution  from  AO  can  be  reduced  or  removed  by 
choosing  some  other  reference  for  measuring  time  differences. 

If  the  arrival  time  at  the  ith  station  is  T^,  referred  to  an 
arbitrary  origin,  then  we  may  fit  a  curve  to  the  data  (as  in 
equation  3)  of  the  form 

T  -  T.  f-  i?  4  +  t  -CD 
u  (8) 

where,  as  before,  b  ^  is  the  distance  from  AO  along  the  azimuth 
to  the  event.  Tq  may  be  determined  from  the  least  squares  fit 
and  used  as  a  reference  time.  The  travel  time  anomalies  with 


SUB  ARRAY  FI 
AZIMUTH  RANGE  290-330 
RESIDUALS  WITH  RESPECT  TO 
EARTH  MODEL  A,  USING 
SUB  ARRAY  AO  AS  TIME  REFERENCE 


(SQN003S)  SIMM  sat) 


DISTANCE  (degrees) 


SUB  ARRAY  FI 


-  42  - 


o 

t^\ 

I 

O 

D 

CNI 


CJ 

UJ 

CO¬ 

CO 

UJ 

oc 


CD 

ct: 


r^j 

<c 


co 

_ i 

<c 

HO 

Q 

CO 

UJ 

oc 


CD  UJ 

z:  cj 
*— <  z; 

CO  UJ 

=>  o; 

LU 
'  LU 
<C  LU 
QJ 


s: 

o  — 


in  co 
i—  <c 
cu 
<3; 


••  •  • •• 


•  • 


•  •  * 


•v.  . 

// 


••• 


••I  • 

1  ,•••«*#•• 


*•  -n 


••• 

•  n  •  . 

•  ••  • 


•  •  •  *  *  #*  • 


r 


• 

•  • 
•  • 


•  ’  ••  -?r 


•  • 


•  • 


«*'#:}  . 


•—  •  •• 

.  1 

•’  •  *•••••  •• 


o 

o 


O 

CD 


O 

OO 


O 

r-\ 


o 

CD 


o 

LT\ 


o 

-zr 


o 


1  1 _ I _ J _ I _ I _ I _ L 

in  o 

o  o 

(SQN003S)  swnais3H 


DISTANCE  (degrees) 


-  43  - 


respect  to  this  reference  for  sub-array  Fl  and  velocity  model  A 
are  shown  in  Figure  24.  There  is  a  considerable  difference 
between  this  graph  and  Figure  23.  Presumably  at  least  part  of 
the  AO  contribution  has  been  removed. 

Both  of  the  techniques  mentioned  above  may  be  criticized. 

The  reality  of  model  A  is  still  questionable,  and  determination 
of  Tp  depends  on  all  of  the  arrivals  at  the  array  and  may  be 
influenced  by  large-scale  structure.  However,  these  techniques 
should  produce  some  improvement  in  the  travel-time  anomalies. 

(ii)  Residuals  from  Curve-Fitting. 

Having  fitted  equation  8  to  the  data,  the  residuals  at.  the 
sub-arrays  may  be  regarded  as  travel-time  anomalies.  This 
definition  of  travel-time  anomalies  has  the  advantage  that  no 
assumptions  about  mantle  structure  are  necessary.  It  has  the 
same  disadvantage  as  before,  that  the  meaning  of  a  fit  such  as 
this  is  not  clear,  and  it  may  be  influenced  by  large  scale 
structure  beneath  the  array. 

Iyer  and  Healy  (1971)  proposed  a  similar  definition  for 
travel  time  anomalies,  though  they  fitted  a  first  degree  equation 
(essentially  equation  2) .  It  turns  out  that  there  are  only 
small  differences  between  the  1st  degree  and  2nd  degree  residuals. 
Iyer  and  Healy  then  averaged  the  residuals  from  9  events  at 
various  azimuths  and  distances,  and  used  the  variation  of  the 
average  residual  across  the  array  to  make  conclusions  about  the 


crustal  structure. 
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Since  these  residuals  have  not  been  studied  in  detail,  the 
2nd  degree  residuals  for  the  9  sub-arrays  used  in  this  study  are 
shown  in  Figures  25-33.  Clearly  there  are  strong  variations  of 
the  residuals  with  epicentral  distance,  at  some  of  the  sub-arrays. 
Similar  complex  variations  occur  as  a  function  of  azimuth.  These 
variations  are  of  the  same  order  of  magnitude  as  those  found  by 
Chiburis  (1966)  ,  but  have  a  different  form  since  they  are  defined 
differently.  The  apparently  small  residuals  at  sub-arrays  F2 
and  F4  are  an  inevitable  result  of  any  fitting  process  to  events 
at  this  azimuth,  since  S  is  largest  for  these  sub-arrays.  The 
magnitude  of  the  variations  throws  some  doubt  on  the  averaging 
method  of  Iyer  and  Healy . 

The  2nd  degree  residual  for  sub-array  FI  (Figure  30)  is 
very  similar  to  the  anomalies  with  respect  to  model  A  (Figures 
24).  This  is  hardly  surprising,  since  model  A  was  determined 
from  the  present  data.  However,  this  points  out  the  underlying 
assumption  of  using  least  square  residuals  as  travel  time 
anomalies.  A  velocity  structure  very  like  model  A  is  implied 
to  be  valid. 

These  residuals  may  be  contoured  acicss  the  array,  in  the 
same  way  that  the  authors  mentioned  have  contoured  travel  time 
anomalies.  Figure  34  shows  an  outline  of  the  array,  and 
Figures  35-38  show  the  mean  residuals  at  distances  of  40,  55, 

70  and  90°,  at  the  same  scale.  The  contours,  though  schematic, 
have  several  interesting  features.  The  values  at' 90°  indicate 
a  Northeast-Southwest  trend  which  is  similar  to  the  conclusions 
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of  others  (e.g.,  Greenfield  and  Sheppard,  1968).  At  closer  dis¬ 
tances  there  appears  to  be  evidence  for  a  superimposed  Northwest- 
Southeast  trend.  The  reality  of  this  trend  and  its  interpretation 
are  not  clear  at  present.  However,  the  variation  in  the  contours 
with  distance  is  significant,  and  has  been  disregarded  in  previous 
studies . 
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Fig.  3A;  ARRAY  GEOMETRY,  AT  SAME 
SCALE  AS  FIGURES  35-38. 


i 


Fig,  37:  CONTOURED  2nd  DEGREE  RESIDUALS  (in 
seconds)  FOR  EVENTS  IN  AZIMUTH  RANGE  290-330, 
AT  AN  EPICENTRAL  DISTANCE  OF  70°. 

FOR  SCALE  SEE  FIGURE  34, 
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Fig,  38:  CONTOURED  2nd  DEGREE  RESIDUALS  (in 
seconds)  FOR  EVENTS  IN  AZIMUTH  RANGE  290-330, 
AT  AN  EP I CENTRAL  DISTANCE  OF  90°, 

FOR  SCALE  SEE  FIGURE  34, 
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VI I .  INTERPRETATION  OF  dT/dA  VERSUS  DISTANCE  (AZIMUTH  RANGE  290-330) 

An  interpretation  of  the  measurements  of  dT/dA  shown  in 
Figures  5  and  6,  for  the  azimuth  range  290-330,  has  been  given  in 
Chinnery  (1'69).  The  inferred  velocity  structures  are  shown  in 
Figures  39-42.  The  velocity  structures  were  obtained  by  a  two- 
stage  process.  First,  a  curve  was  fitted  visually  to  the  dT/dA 
data,  and  then  this  was  inverted  into  a  velocity  structure  using 
the  well-known  Weichert-Herglotz  technique.’  The  resulting 
velocity  structure  was  then  used  to  recalculate  dT/dA  as  a 
function  of  distance,  and  it  is  this  curve  that  is  shown  through 
the  data  points.  This  technique  largely  overcomes  any  uniqueneso 
problems  in  the  Wei cKert-Herg lots  method,  except  those  introduced 
by  scatter  in  the  data.  The  additional  data  points  shown  in 
Figures  5  and  6  have  not  made  any  significant  change  to  this 
interpretation. 

There  are  several  points  about  this  interpretation  that 
deserve  consideration  here,  in  the  light  of  data  obtained  since 
the  paper  was  published.  First,  some  doubts  were  thrown  on  the 
interpretation  because  station  corrections  at  the  array  were  not 
included.  The  reason  for  this  given  in  the  paper  was  that  the 
addition  of  station  corrections  will  modify  the  shape  of  the 
curve  of  dT/dA  against  distance,  and  so  change  the  travel  time 
residuals  of  the  curve  with  respect  to  the  Jef freys-Bullen  table. 
(These  residuals  are  simply  the  cumulative  area  between  the  two 
curves.)  The  nature  of  the  modification  will  depend  on  the 
structure  assumed  in  the  calculation  of  the  station  corrections. 


dT/dA  (i*o/d*g) 


dT/dA  (tecs /deg) 


A  simple  dipping  interface  will  raise  or  lower  the  observations 
with  respect  to  the  Jef freys-Bullen  curve.  More  complicated 
structures  will  distort  the  curve.  However,  the  T-JB  residuals 
of  the  uncorrected  observations  agree  surprisingly  well  with  the 
T-JB  residuals  at  North  American  stations  from  the  Longshot 
nuclear  event,  in  the  Aleutian  Islands  (Figure  43) .  The  travel 
paths  for  this  event  were  very  similar  to  those  involved  in  this 
study.  It  is  therefore  concluded  that  station  corrections  for 
the  azimuth  range  290-330°  must  be  very  small,  and  may  to  a  first 
approximation  be  neglected. 

More  recently  a  number  of  other  arrays  have  been  used  to 
determine  dT/d  A  •  Figure  44  shows  that  results  of  Johnson 
(1969)  ,  who  used  the  TFSO  array,  and  Corbishley  (1970) ,  who 
combined  the  data  from  four  arrays.  Figure  45  shows  the  results 
obtained  by  Wright  (1971)  using  the  Warramunga  array  in  Australia. 
Although  none  of  these  newer  studies  have  the  resolution  (small 
scatter)  and  amount  of  data  contained  in  Figures  5  and  6 ,  the 
similarities  between  the  results  are  striking.  Again,  this 
strongly  argues  that  the  data  in  the  present  study  should  not  be 
distorted  by  the  addition  of  station  corrections. 

Further  evidence  on  this  point  is  obtained  from  amplitude 

studies,  since  the  amplitude  of  an  srrival  det ?nds  strongly  on 

2  i 

the  second  derivative  d  T/d A  .  Figure  46  shows  how  the  logarithm 
of  the  slope  of  the  dT/d A  graph  varies  with  distance,  compared 
with  the  amplitude  observations  of  Carpenter  et^  al_.  (1967),  using 
an  arbitrary  zero.  Again  a  remarkable  agreement  is  seen. 
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Fig.  43:  Travel  time  residuals  with  respect  tc  the  Jeffreys 
Sullen  Tables;  above,  observations  of  the  LONGSHOT  test, 
below,  corresponding  to  models  A  and  B  in  figures  39-42. 
(After  Chinnery,  1969) 


-  67  - 


Fig.  44:  Other  observations  of  dT/dA  ;  above,  from  Johnson 
(1969),  below,  from  Corbishley  (1970). 
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Fig.  45:  Observations  of  dT/dA  by  ’/right  (1971) 
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We  conclude  that  the  smooth  model  A  or  the  less  smooth 


model  B,  shown  in  Figures  39-42,  are  reliable  structures  for 
the  paths  concerned  in  this  study. 

We  may  also  return  to  our  earlier  comment  that  the  scatter 
of  the  dT/dA  observations  is  unusually  large  in  two  distance 
ranges,  32-37°  and  83-88°.  The  latter  will  be  discussed  in  the 
next  section.  The  former  has  become  interesting  because  of  a 
suggestion,  by  Wright  (1968,  1971),  that  a  triplication  may  occur 
in  the  vicinity  of  35°,  as  the  result  of  a  low  velocity  layer  in 
the  mantle  at  a  depth  of  about  850  km.  Wrights  dT/dA  data  were 
obtained  using  the  Warramunga  array,  and  show  a  very  similar 
structure  to  those  in  Figure  5  at  this  distance.  It  is  very 
likely,  therefore,  that  this  effect  is  real.  A  possible  inter¬ 
pretation  of  the  present  data,  in  terms  of  a  low  velocity  layer, 
is  shown  in  Figure  47.  A  variety  of  models  are  possible,  since 
data  in  the  vicinity  of  35°  is  sparse,  but  the  model  AN  shown 
is  representative.  Further  studies  of  this  feature  are  continuing. 
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VIII.  AZIMUTHAL  VARIATIONS 


Having  established  an  apparently  self-consistent  interpretation 
for  the  data  from  the  azimuth  range  290-330,  it  is  natural  to 
attempt  the  same  process  for  data  from  other  azimuths.  It  is  at 
this  point  that  the  differences  between  the  observations  at  dif¬ 
ferent  azimuths  become  most  apparent. 

Fairly  complete  information  as  a  function  of  distance  is 
available  for  events  in  the  azimuth  range  130-170  (Figures  11  and 
12) .  In  order  to  make  a  direct  comparison  with  the  290-330  data, 
a  visual  fit  to  the  data  in  Figures  11  and  12  is  shown,  together 
with  model  A  (see  previous  section)  and  the  Jef freys-Bullen 
curve,  in  Figure  48.  The  differences  between  the  two  sets  of 
data  are  extreme.  This  may  be  a  common  feature  of  array  measure¬ 
ments;  certainly  similar  differences  were  observed  at  the 
Warramunga  array  by  Wright  (1971),  as  is  shown  in  Figure  49. 

There  are  three  possible  ways  in  which  these  differences  may 
occur;  they  may  be  due  to  crustal  effects  at  the  array,  or  to 
lateral  variations  in  mantle  structure,  or  to  effects  in  the 
vicinity  of  the  source  (or,  most  likely,  to  some  combination 
of  all  three) .  These  possibilities  are  discussed  separately 
below. 

(i)  Effects  of  Crustal  Structure. 

Crustal  structure  at  the  array  must  be  a  major  contributor 
to  the  differences  between  the  data  from  these  two  azimuth 
ranges.  However,  the  separation  of  this  contribution  from  the 
other  two  possibilities  is  not  easy. 


dT/dA.Mc/dag 
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Fig.  49:  Comparison  of  the  data  observed  in  two  directions 
from  the  V/arramunga  array  (after  './right,  1971). 


-  75  - 


If  the  paths  from  130-170  traversed  a  simple  dipping  structure 
in  the  vicinity  of  the  array,  the  dT/dA  curve  would  be  displaced 
parallel  to  the  dT/dA  axis.  It  is  interesting  to  note  that  the 
curves  from  the  two  azimuth  ranges  are  quite  parallel  in  the 

O  0 

distance  range  30-55  ,  with  the  130-170  data  displaced  upward  by 
0.5  seconds/degree.  This  is  a  substantial  offset,  and  indicates 
a  comparatively  steeping  dipping  interface  with  considerable 
velocity  contrast  (Niazi,  1966). 

The  parallelism  disappears  at  distances  greater  than  about 
55°,  and  it  becomes  much  more  difficult  (though  presumably  not 
impossible)  to  relate  the  differences  to  a  specific  crustal 
model. 

(ii)  Lateral  Variations  in  Mantle  Structure. 

The  problem  is  compounded  when  these  curves  are  studied  in 

o 

detail.  There  are  a  number  of  anomalous  regions  in  the  290-330 
curve  that  appear  to  be  due  to  inhomogeneities  in  the  mantle 
(Toksoz  et  al . ,  1967).  None  of  these,  including  the  feature 

O 

shown  in  Figure  47/  are  apparent  in  the  130-170  data.  This  is 

supported  by  the  d  T/dA  observations  given  in  section  4  of 

this  report.  On  the  other  hand,  there  is  a  variety  of  evidence 

o 

to  suggest  that  the  interpretation  of  the  290-330  data  is 
reasonable  (see  preceding  section) .  We  are  therefore  reluctantly 
led  to  the  conclusion  that  the  seismic  paths  from  the  South 
American  earthquakes  traverse  a  highly  anomalous  section  of  the 
mantle  (Chinnery,  1967). 
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The  degree  of  anomaly  implied  by  these  data  is  very  large. 

o 

The  flattened  portion  of  the  290-330  dT/dA  curve  at  about  70° 

suggests  the  presence  of  a  phase  transition  at  a  depth  of  about 

2000  km  in  the  mantle  (Chinnery,  1969b).  The  continuity  of  the 

130-170  data  in  this  region  appears  to  preclude  the  presence 

of  such  a  feature  anywhere  in  the  depth  range  1500-2500  km.  On 

the  other  hand,  the  data  from  the  vicinity  of  Easter  Island 

u 

(azimuth  range  170-200,  Figure  14)  do  appear  to  indicate  the 
presence  of  this  feature . 

It  is  not  clear  whether  lateral  variations  in  deep  mantle 

structure  of  this  magnitude  are  permissible.  Certain  types  of 

convective  flow  in  the  mantle  might  lead  to  variations  like 

this,  but  the  present  data  have  sufficient  ambiguity  that  they 

cannot  prove  the  existence  of  such  features.  It  is  even  possible 

that  the  presence  of  the  phase  transitions  is  obscured  in  the 
© 

130-170  data  by  a  complex  array  structure,  though  this  seems 
unlikely . 

One  is  led  to  the  following  conclusions.  If  the  inter- 

o 

pretation  of  the  290-330  data  is  valid,  then  seismic  rays  from 

this  direction  see  a  very  simple  structure  at  the  array  (on 

o 

the  average).  Then,  the  130-170  data  imply  either  an  enormously 
complicated  crustal  structure  at  the  array  for  paths  from  this 
azimuth,  or  major  lateral  variations  in  mantle  velocity  structure. 
Neither  of  the  possibilities  are  very  desirable.  The  only 
alternative  is  that  the  crustal  structure  at  the  array  is  so 
complex  that  a  "noise"  is  added  to  the  observations  at  all 
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azimuths.  In  this  case  it  is  surprising  that  the  interpretation 
0 

of  the  290-330  data  agrees  so  well  with  ether  studies. 

(iii)  Source  Structure. 

The  possibility  of  contamination  of  dT/dA  data  by  structure 
close  to  the  source  has  received  little  attention  in  the  past. 

This  is  not  surprising,  since  the  angular  separation  of  the  rays 
that  will  ultimately  arrive  at  the  extremes  of  the  array  is  very 
small  at  the  source  (much  less  than  1°  for  any  mantle  structure) . 
Only  an  extraordinary  structure  could  introduce  a  significant 
time  difference  into  two  rays  so  close  together.  However,  there 
is  at  least  one  instance  in  the  present  data  where  this  appears 
to  have  occurred. 

Figure  50  shows  an  enlarged  view  of  the  scattered  data  at 

o  ® 

distances  of  80-90  ,  in  the  azimuth  range  290-330  (Figure  6). 

The  data  have  been  studied  carefully,  and  all  doubtful  information 
removed.  At  about  80°  from  the  array,  the  distribution  of  earth¬ 
quakes  bifurcates  at  the  well-known  junction  just  south  of  Japan. 
The  data  marked  by  solid  circles  and  triangles  are  from  events 
in  Southern  Honshu,  Kyushu,  and  the  Ryukyu  Islands.  Those  marked 
by  open  circles  and  triangles  are  from  events  in  the  Bonin  and 
Mariana  Islands.  The  dT/dA  values  from  these  two  regions,  at 
the  same  distance,  differ  by  as  much  as  0.3  seconds/degree. 

Since  the  azimuths  concerned  are  so  similar,  it  is  unlikely 
that  this  could  be  produced  by  array  structure,  and  this  is 
substantiated  by  the  agreement  between  the  sets  of  data  at 
distances  less  than  80°  and  greater  than  90°. 
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At  first  sight,  it  is  likely  that  this  difference  is  the  result 
of  source  structure.  The  azimuth  to  the  array  from  events  in  the 
Ryukyu  Islands  is  parallel  to  the  oceanic  trench.  However,  the 
azimuth  to  the  array  lies  at  a  large  angle  to  the  island  arc 
marked  by  the  Bonin  and  Mariana  Islands.  However,  on  closer 
study,  it  is  not  easy  to  account  for  the  observed  differences  in 
terms  of  normal  island  arc  structures,  for  the  dip  of  the  Benioff 
zone  is  away  from  the  direction  toward  the  array.  This  means  that 
the  difference  must  arise  either  in  the  lithosphere  (all  the  events 
are  shallow,  focal  depths  less  than  100  km) ,  or  in  the  upper  mantle 
on  the  ocean  side  of  the  downgoing  slab.  Since  focal  depth  appears 
to  have  no  significant  influence  on  the  measured  value  of  dT/dA  , 
it  is  not  clear  that  the  former  possibility  is  reasonable.  However, 
if  mantle  structure  is  the  cause,  large  lateral  variations  in 
velocity  are  implied.  Until  more  is  known  about  the  structure  of 
the  mantle  in  the  vicinity  of  a  downgoing  slab,  it  does  not  seem 
possible  to  produce  a  consistent  interpretation  of  these  data. 

It  appears,  then,  that  array  structure,  lateral  variations 
in  mantle  structure,  and  source  structure,  may  each  at  times 
contribute  significantly  to  contamination  of  dT/dA  measurements. 
There  are  severe  problems  in  attempting  to  separate  out  anyone 
of  the  three  for  individual  study. 
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CONCLUDING  REMARKS 

The  data  set  that  has  resulted  from  this  study  consists  of 
the  following: 

(i)  P-wave  arrival  time  information  at  the  9  subarrays 
described  for  the  more  than  3000  events  studied. 

These  are  punched  on  one  card  for  each  event,  and 
included  are  an  event  code,  and  the  USCGS  latitude, 
longitude  and  depth  of  the  event. 

(ii)  A  one  page  computer  print-out,  as  shown  in  Figure  2b, 
for  each  event. 

The  following  programs  have  been  written  in  Fortran  IV 
for  the  IBM  360  computer: 

(i)  EVENTS:  a  versatile  program  for  the  analysis  of  the 
travel  time  data,  giving  the  output  shown  in  Figure  2b. 

(ii)  TIMDIS :  a  program  that  will  accept  tiny  velocity-depth 
model,  interpolate  between  data  points  to  improve 
continuity,  and  then  calculate,  for  a  given  value  of 
dT/dA,  a  variety  of  travel  time  information.  This 
includes  travel  time,  distance,  depth  of  lowest  point 
of  ray  path,  residual  from  the  Jef f reys-Bullen  travel 
time  table,  the  second  derivative  of  the  travel  time 
curve,  PcP  times,  etc. 

(iii)  WIEHER:  a  program  for  the  inversion  of  dT/d&  data  to 
obtain  a  velocity-depth  model,  using  the  Wiechert- 
Herglotz  method. 

The  graphs  shown  in  this  report  are  only  a  small  fraction  of 


those  that  could  be  drawn  using  the  existing  information. 
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FINAL  REPORT 

Contract  F44620-68-C-0082 

The  results  of  research  performed  under  this  contract  fall 
into  the  following  categories: 

1.  Tf.sk  A:  Investigations  of  P-wave  Arrivals  at  LASA. 

This  research  is  discussed  fully  in  the  attached 
technical  report. 

2 .  Task  B:  Research  in  Acoustic-Gravity  !,Taves. 

This  research  was  carried  out  by  Professor  David  G. 
Harkrider  (now  at  California  Institute  of  Technology) .  A 
final  technical  report  of  this  research  has  been  submitted 
separate ly . 

3 .  Research  During  the  Period  8/1/70  to  7/31/71. 

A  no-cost  extension  of  the  contract  was  granted  for 
the  period  8/1/70  to  7/31/71,  under  Amendment/Modification 
No:  P002.  Research  carried  out  during  this  period  is 
described  below.  Much  of  this  research  is  still  in  progress. 
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Research  Under  Amendment /Hod if ication  No.  P002 

This  research  has  been  concentrated  on  "Continuing  Research 
Objective  M"  (see  original  proposal,  page  9),  which  is  concerned 
with  the  residual  displacements  that  accompany  earthquake  sources. 
Portions  of  this  researcn  are  described  in  three  technical  papers 
(one  in  press,  and  two  in  the  late  stages  of  preparation). 

Abstracts  of  these  papers  are  given  below: 

1.  On  the  Correlation  between  Earthquake  Occurrence  and  Dis¬ 
turbances  in  the  Path  of  the  Rotation  Pole,  by  M.  A.  Chinnery 
and  F.  J.  Wells.  (Given  orally  at  the  International  Astro¬ 
nomical  Union  Symposium  #48  on  Rotation  of  the  Earth, 

Morioka,  Japan,  May,  1971;  paper  will  appear  in  the  pro¬ 
ceedings  of  the  symposium.) 

The  hypothesis  that  earthquakes  may  be  the  principal 
-'xcitation  of  the  Chandler  motion  of  the  rotation  pole  is 
examined  in  the  light  of  recent  theoretical  and  observational 
developments.  There  is  some  doubt  about  the  amount  of 
excitation  by  a  large  earthquake  necessary  to  maintain  the 
Chandler  Wobble,  but  it  appears  to  be  about  10  feet.  Theo¬ 
retical  calculations  for  the  Alaskan  Earthquake  (M  =  8%) 
give  available  excitations  in  the  range  1-5  feet,  but  there 
are  considerable  uncertainties  in  these  calculations.  Earth¬ 
quakes  may  be  able  to  provide  all  of  the  required  excitation, 
or  only  a  small  portion  (10%  or  less) .  The  problem  is  con¬ 
fused  by  observational  studies,  which  show  differences  between 
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various  sets  of  data  on  polar  motion  which  seem  to  be  larger 
than  the  expected  error  in  each  set.  The  earthquake  hypo¬ 
thesis,  though  reasonable,  is  still  very  much  open  to  debate. 

2.  The  Inertia  Changes  due  to  Faulting  in  a  Homogeneous  Sphere, 
by  M.  A.  Chinnery  and  J.  R.  Pice  (manuscript,  in  late  stages 
of  preparation) . 

A  very  simple  and  elegant  theory  for  the  changes  in 
the  inertia  tensor  due  to  slip  over  an  arbitrary  surface 
inside  a  homogeneous  sphere  is  presented.  Full  expressions 
are  given  for  the  case  of  a  point  slip  source,  for  the 
inertia  changes  and  AI33*  These  results  are 

compared  with  those  of  Ben-Menahem  and  co-workers,  and 
Dahlen. 

3.  Two-Dimensional  Faulting  in  a  Layered  Half-space,  by 

M.  A.  Chinnery  and  D.  B.  Jovanovich.  (Presented  orally  at 
the  April,  1971  meetings  of  the  American  Geophysical  Union, 
Washington,  D.  C.;  manuscript  in  late  stages  of  preparation). 

The  method  of  images  is  used  to  evaluate  the  displace¬ 
ment  field  due  to  a  very  long  (2-dimensional)  fault  in  two 
layers  over  a  half  space,  with  arbitrary  rigidities.  The 
model  is  used  to  investigate  the  effect  on  earthquake 
displacements  of  the  lithosphere-asthenosphere  boundary, 
the  low  velocity  layer,  and  possible  thin  soft  layers  in 
the  upper  mantle.  Each  of  these  produce  an  amplification 
of  fhc  curfarc  displacements  over  a  uniform  half-space 
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model.  It  is  shown  that  thin  layers  of  soft  material  may 
effectively  decouple  earthquake  displacements  from  the 
earth's  interior. 


Continuing  Research 

Continuing  research,  which  will  be  reported  later,  and 
which  will  acknowledge  the  support  received  under  this  contract, 
is  described  below. 

Considerations  such  as  those  contained  in  reference  #3 
above  have  strengthened  our  interest  in  the  far  displacement 
fields  of  earthquakes.  It  is  likely  that  these  fields  may 
contain  significant  information  about  source  mechanism  and  earth 
structure,  including  some  information  not  readily  available  from 
seismic  wave  studies. 

It  is  impossible,  at  present,  to  detect  these  fields 
directly.  Strain  step  observations  permit  an  indirect  method 
of  detection,  but  relatively  few  strain  meters  are  available, 
and  the  observations  of  strain  steps  are  contaminated  to  an 
unknown  extent  by  local  structure  and  instrumental  effects. 

Astronomical  observations  of  latitude  variation  are  one 
possible  way  to  detect  the  total  mass  redistribution  associated 
with  these  displacement  fields.  In  particular,  it  has  been 
suggested  that 

(i)  The  total  displacement  field  of  a  large  earthquake 
may  be  sufficient  to  cause  a  reasonable  change  in  the  moments  of 
inertia  of  the  earth,  and  so  excite  the  Chandler  Wobble. 


(ii)  This  excitation  may  precede  the  occurrence  of  the 
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earthquake  by  several  days. 

(iii)  A  natural  corollary  to  these  is  that  information 
about  large-scale  motions  in  the  earth  (apart  from  earthquakes) 
may  be  obtained  from  a  detailed  study  of  the  path  of  the  rotation 
pole . 

We  have  therefore  been  led  to  examine,  in  some  detail, 
the  astronomical  observations  of  the  path  of  the  rotation  pole. 

It  turns  out  that  these  observations  are  extremely  complex,  nnd 
are  probably  much  less  reliable  than  has  been  reported  in  the 
literature. 

The  most  serious  problem  appears  to  be  the  separation  of 
the  free  Eulerian  motion  of  the  earth  (14  month  period)  from  a 
variety  of  effects  with  a  12  month  period.  The  latter  include 
the  forced  motion  of  the  earth  due  to  atmospheric  motions, 
together  with  a  series  of  local  effects  (local  meteorology, 
star  catalogue  errors,  operator  changes,  etc.)  that  each  con¬ 
tribute  substantially  (and  with  variable  phase)  to  the  12  month 
component.  The  separation  of  these  two  terms  is  essential  to 
any  study  of  the  excitation  of  the  free  motion  b”  earthquakes. 
However,  examination  of  the  published  techniques  used  to  carry 
out  the  separation  has  shown  that  they  are  quite  inadequate.. 

This  arises  partly  because  of  the  small  separation  between  the 
two  spectral  peaks,  and  partly  because  of  the  variability  of 
amplitude  and  phase  of  the  12  month  component. 

Present  investigations  are  concentrating  on  the  statistical 
characteristics  of  these  data,  and  the  extent  to  which  separation 
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is  possible.  This  in  turn  will  lead  to  a  more  complete  under¬ 
standing  of  the  validity  of  the  hypotheses  mentioned  above.  It 
is  not  expected  that  definitive  answers  to  these  questions  can 
be  obtained  from  the  presently  available  data.  However,  within 
a  few  years  it  is  expected  that  satellite  data  will  become 
available  that  will  permit  much  more  accurate  (by  an  order  of 
magnitude)  determinations  of  pole  position.  This  study  is 
expected  to  lay  the  groundwork  for  the  analysis  of  this  new  data. 


